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ABSTRACT 



This thesis investigates the effects of atmospheric 
turbulence on the spatial properties of the phase of transmitted 
beams. Experimental measurements in the form of a series of 
photographs of interference patterns between two beams were 
made using a modified optical heterodyne system. The two beams 
of the heterodyne system were mixed with a very small, control- 
lable angular difference between their directions of approach to 
the detector. The resulting fringe patterns were effected by the 
atmosphere as differences between succeeding photographs show. 

The first order effect of the atmosphere was found to be a tilt 
relative to the direction of travel. Tilt differences as small as 
5 seconds of arc can be measured using this system with assurance; 
smaller values could be confused with systematic effects. 
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INTRODUCTION 



The advent of the laser has stimulated interest in communl- 
cation at ootical frequencies. The potential bandwidths and an- 
tenna gains hold great promise, but the degrading effects of the 
atmosphere on propagating optical wavelength beams have prevented 
immediate widespread usage. Present efforts in the area of atmos- 
pheric optical communication are mainly directed towards gaining 
an understanding of the atmospheric phenomena and characterizing 
the received fields after propagation through the atmosphere. 

This information is necessary to build a statistical model of the 
channel at optical and near infrared frequencies. A statistical 
model will allow the use of tools from modern communication theory 
to establish the ultimate limits of optical communication links 
utilizing atmospheric channels. 

Atmospheric effects on propagated fields at optical and near 

infrared frequencies can be broadly put into ttree catagorles! 

1 . Molecular absorption 

2 . Particle scattering 

3 . Turbulence 

The first two are more or less suggestive and under some condi- 
tions these two effects can be lumped as a path-length dependent 
loss term. This assumption depends strongly on atmospheric condi- 
tions and propagating wavelengths. However, with a clear atmosphere 
this model is reasonable at a wavelength of 6328A. This thesis 
deals with the third effect, turbulence. 

Turbulence is the name used to describe the effects of small 
variations In the local refractive index caused by turbulent mix- 
ing of the atmosphere under wind and temperature gradients. At- 
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rnospheric turbulence results In an inhomogeneous medium, and pro- 
pagation through this medium causes random spatial and temporal 
fluctuations in amplitude and phase of the received field. A 
body of work has been done on the amplitude fluctuations, called 
scintillation, with beams. There has been less emphasis on the 
phase process due to the difficulty of phase measurements. 

This thesis considers a method of phase measurement using 
an optical heterodyne system and evaluation of the experimental 
measurements. 

Chapter 1 presents a brief survey of the pertinent theory 
and modeling necessary for this thesis. 

Chapter 2 describes the experimental system used to hetero- 
dyne at 6328 ^ and the data analysis system in general. Appendix 
A contains a more detailed explanation of the system, special 
devices, and electronic circuits. 

^Chapter 3 presents the experimental results. Appendix B con 
tains more detailed information on the method of data extraction. 

Chapter 4 considers the experimental results and conclusions 



CHAPTER 1 

This chapter is concerned with some of the theoretical re- 
suits developed for beam propagation through the turbulent atmos- 
where with emphasis on those results pertinent to the experimental 
system constructed. 

1-1 Atmospheric Turbulence 

The effects of atmospheric turbulence on optical wavelengths 
have long been known to astronomers; the twinkling of starlight on 
a windy night is a common example. The possible use of the laser 
for atmospheric communication links has renewed interest in the 

| 

study of atmospheric turbulence and its effects. 

Turbulence is small scale vatiations of the index of re- 
fraction, typically on the order of 1 part in 10^, caused by 
temperature discontinuities set up by turbulent mixing of the at- 
mosphere. These refractive index changes were first described using 
a Gaussian correlation function which was an analytically convenient 
method best led to questionable results. A second path taken by 
Tatarski applied the Kolmogorov theory of turbulent flow at very 
large Reynolds numbers to the index variations^ This analysis 
leads to a second moment description of the medium, a method that 
is not alone sufficient to describe atmospheric turbulence. How- 
ever, the fields propagated, through the inhomogeneous medium will 
have lognormal statistics; thus a second moment description of the 
medium will be adequate to handle the received field statistics. 

The results of Tatarskl's work show that the turbulence is locally 
isotropic and homogeneous and has a structure function of the form 1 
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2/3 



f- o 4p£Lo 



D n (p,R)= < 




( 1 ) 




where 



D n (p,R)= [n(r. ± ) - n( r-^+p ) ] 2 



( 2 ) 




where R Is a vector. A bar will be 
used to denote vectors. 



P 

2 



P 




constant associated with the magnitude of the 
turbulence at R 



n(r) = refractive index at r 

^ o = inner scale of turbulence; on the order of 1mm 
and related to the breakup of_JLaminar flow. At 
distances greater than^o, n(r) tends to be dis- 
continuous » 

Lo = outer scale of turbulence; on the order of meters 
and related to the gross layer structure of the 
atmosphere or the distance between obstacles. 



The models while not as convenient analytically, is generally 

description 

considered better than the Gaussian correlation functior/and has 



1-2 Pro paga tion thr o u gh the^ t urbulen t atmosphere 

The available theoretical solutions generally are for Infinite 
plane wave, sources and assume a "frozen" atmosphere, where the 
turbulence is fixed in space. The frozen-in structure is assumed 
to move across the transmission path due to the component of wind 
normal to the propagation path to account for the temporal behavior. 

The three basic approaches to the turbulent propagation problem arel 

1 

using a Rytov approximation to the wave equation, solving the wave 



been the basis for most further work 
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2 

equation for the mutual coherence function, and a geometric optics 

3 

type of solution. The last method, geometrical optics, is general- 
ly considered inferior since it does not account for diffraction 
effects on long paths with small scale inhomogenetics. 

1-3 Received Field Parameters 



Using the '’frozen" atmosphere assumption, the central limit 

theorem, and one of the solution approaches, it can be shown that 

a time sample of the received field envelope is a lognormal random 

variable, and that a time sample of the received field phase is a 

Gaussian random variable. There is evidence to expect that the 

form of these statistics is unchanged for beam propagation. 

Another measure of importance is the spatial coherence of the 

phase of the received field. Tatarski has derived an expression for 

the case of a linearly polarized plane wave with the atmospheric 

2 _ 

turbulence described as in Eq. 1, and C (R) a constant independent 
_ 1 n 
of R . Then 

225/3 

Ity(p) - aC n k Lp (3) 

where D^>(p) is the phase structure function, k the wave number of 
the radiation, L the path length, and 



:.91 



P > XL 



a = 

1.45 ?o4< p« XL 

Using the expression for D ( £(p’) in Eq. 3, a coherence diameter for 

some selected, mean-square difference in phase can be calculated. 

9 4 

Qne example quoted in Hoversten from Goldstein is 



D 



coherent 



= 1.75 x 10" 



°n L 



(4) 
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Fried has used a form of Eq. 3 to show that wave^front tilting is 

5 

the largest contributor to distortion. 

/mother parameter is the frequency spectrum of the amplitude 
and phase turbulence processes. Measurements of Intensity spec- 
truras have shown them to be no more than several hundred cycles 

V 

wide for most conditions. Measurements for the phase process are 

not as plentiful, but recent measurements have shovrn that the spec- 

s 

trum is no more than several hundred, cycles wide , 

1 - Communication Mod el 



A suggested model describing the effects of the turbulent 
atmosphere as a communication channel ,when particle scattering 

is not important, is a multiplicative channel with additive back- 

? 

ground noise. A linearly polarized field (chosen for simplicity) 
in the absence of losses, turbulence and background noise would 
be received as 

<- _ 

E rec (r,t) = Real|E(r,t)e $ (5) 

When turbulence and background noise are present the received field 
(again in the absence of losses) would be 

_ r 3<«t $ ( r , t ) ^n(r,t) Jwt \ 

E rec ( r J t ) = Real | E(r,t)e e + E n (r,t)e e j 

( 6 ) 

where 

e^ ' * ' represents the effect of the turbulent channel 
tf(r,t) = X(r,t) j^(r,t) 

X(r,t) represents the scintillation process 
<£>(r,t) represents the phase distortion process 
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J© n (r,t) 

E (r,t)e represents the relevant polarization 

component of the background radiation; 
a complex Gaussian process with a very 
wide band 'width on both E n and© n . 

1 - 5 Optical Heterodyning 

Two waves simultaneously incident on a surface will cause an 
interference pattern. This can be modeled simply, using plane waves 
of the same frequency and linear polarization. When the two waves 
beat together from slightly different directions, 0^ and a 

fringe pattern whose period is related to the difference in approach 
angle is created. If one of the waves has been sent through the 
atmosphere, it xfill have a phase term proportional to the phase dis- 
tortion process. Thus if complex notation is used, suppressing 
the time variations for simplicity, the fields can be des- 

cribed as (see Fig, 1) 



iW 1 -.*) = El® 



-)Bz X( r,t) 3<kr.t) 



where 



t -JB z-JB x 

E Lo (r,t) = S 2 e z x 



I B 1 1 = l B 2 I =4? + B 2 



(8) 



( 9 ) 



= j~- and 0-j and g ? are 
wave Vectors; 



z ^x 

Then on the surface (z=o) shown in Fig. 1 the total field can be 
written as 



£total (rst) = S rec (^t) + E Lo (r,t) 

X(pi*t) J^(p,t) -jB x 

= F^e e + E 2 e x 



( 10 ) 



where 
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e = '■ x x + V 




Then the intensity on the surface is 
^ [StotalS?otal} 

: , 2 2X(p,t) 2 X(p,t) 

[l-J e + I E| + 2IE.J Egl e cos[B x x-<p( p, t)-cx 

( 11 ) 

where is the impedance of the atmosphere and a represents the 
phase angle between the complex amplitudes and Eg. Thus a por- 
tion of the intensity is a term phase modulated by the atmospheric 
phase distortion. In the experimental section of this thesis a 
system to exploit this effect is described in detail. 




CHAPTER 2 
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This chapter describes the system developed, to create a 
spatial enterference pattern between two beams of coherent 'light, 
one modulated by atmospheric turbulence, 

2- 1 S ystem Layout a n d Equipmen t 

The basic arrangement for heterodyning is shown in Fig. 2. 

The source was a Spectra-Physics Model 119 Helium-Neon gas laser- 
operating at a wavelength of 6328A. The Model 119 laser has a 
rated continuous output power of 0.1 milliwatt although measure- 
ments indicated slightly higher values during experimentation. The 
beam from one end was collimated using a Spectra-Physics Model 331- 
333 spatial filter-collimator and transmitted over an atmospheric 
path. The down range reflector, shown in Fig. 3 and Fig. 4, is a 
sixteen-inch diameter flat mirror mounted in an aluminum collar and 
frame. The collar and frame have precision adjustments through screw- 
arm systems that allow very accurate pointing of the return beam. 

The beam from the opposite end of the laser was reflected off two 
flat first surface mirrors mounted on precision adjustment holders 
to reach the combining beam splitter. No lens vras used in this non- 
atmospheric path to minimize extraneous lens aberrations. As this 
beam vras close to plane wave conditions by the time it reached the 
combining beam splitter, its curvature was easily matched by the 
telescope output using the telescope eyepiece. 

The receiving telescope is a three-inch reflector with an 
eighteen-inch focal length. The lens from a Praktlca FX-3 camera 
(f/2.8, 5° millimeter focal length) vras used as the eyepiece for 
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Figure 2. System Layout for Optical Heterodyne 




Figure 3* Front View of Downrange Reflector and Mount 




Figure 4. Rear View of Downrange Reflector and Mount 
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several reasons; the 50 millimeter focal length meant the telescope 
output beam size closely matched the local oscillator beam size 
and the focusing ring of the eyepiece allowed adjustments to be 
made with a minimum of mechanical movement. 

The two signals were combined using a fifteen-millimeter 
prismatic bea.m splitter mounted on a 1 ine adjustment holder. This 
mount was used to manipulate the fringe spacing and orientation 
since relatively small changes in angular separation can be made 
most easily at this point. 

All of this equipment, with the exception of the large down- 
range reflector, was bolted down on a two—incn thick aluminum op- 
tical bench housed in a penthouse on the roof of Building 20 at M.-I.T. 
(See Figs. 5 and 6). The downrange reflector was housed in a 
protective box on the other end of the roof of Building 20, mak- 
ing the atmospheric path approximately 200 meters rouno. trip. 

2-2 Pat te rn Detecti on 

The Interference pattern of the two beams was recorded on 
black and white film using a Praktica FX~3 single lens reflex 

camera. The 35mm negatives were developed using Kodak formula D-76 
developer, standard for the Tri-X film used, and were mounted in 
holders to be used as slides for later projection. The negatives 
were used directly for this projection step since the pattern of 
light and dark fringes was the same whether or not the image inten- 
sities had been reversed. 

The fringe patterns 'were examined in detail using a slide 
projector to put the image on a screen. This image was then 
scanned using a Packard Bell Model 920 television camera. The 




Figure 5« Optical Bench Layout Showing Laser/Collimator, 

Receiving Telescope and Eyepiece, Combining Beam 
Splitter, and Camera 




Figure 6. Optical Bench Showing Local Oscillator Path 
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video output of the television camera was displayed on a Tech- 
tronix Model 5^5 oscilloscope as described in section 2-*!-. 

2-3 Experimental Procedure 

The alignment between the source and receiver was acpomplishod 
by finding the laser (not turned on) image in the receiving tele- 
scope and down-range mirror system during daylight hours. Thus 
the optical axis between laser, mirror, and telescope wa.s estab- 
lished. The laser beam could by angular adjustment of the laser 
mount then be put on the mirror at night to have this beam direction 
coincide with the optical axis. 

At this point the beam emerging from the opposite end of the 
laser was positioned with the small flat mirrors so that it pointed 
at the center of the beam splitter. All further adjustments were 
made using the beam splitter. The local oscillator signal, making 
a right angle change in direction was much more sensitive to the 
movement of the beam splitter than the beam coming out of the tele- 
scope going straight through. Thus the local oscillator beam could 
be moved without disturbing the telescope output. 

The beams coming out of the beam splitter could be safely ob- 
served by eye after they had passed through a Kodak# 2 Wratten 
filter. Without an eyepiece the beams appeared as pencil points 
which could be adjusted to coincide. An eyepiece could then be 
used to see the resulting fringe pattern. 

After establishing a fringe pattern, the filter was removed, 
and the Praktica FX~3 camera without lens was used to record se- 
quences of three to six consecutive frames of the same fringe pat- 
tern ie., three to six frames were exposed without altering the 
optical system. 



2 - 4 Data Presentati on 
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I The slides used for analysis were scanned using a Packard 

i 

Bell 920 television camera synchronised to an external clock. 

'This same clock output v.’as used to trigger the oscilloscope. 

The video output was displayed in two distinctive manners. 

The first method was designed to obtain an overall picture of the 
fringe process across the aperture , the second method to look at 
narrow strips of the fringe process across the aperture. 

The first method is shown in block diagram form in Figure 7* 
■The horizontal sweep lines, modulated by the alternating dark and 
light strips of the projected fringe patterns, were displayed one 
on top of the other by setting the time base of the oscilloscope 
to the horizontal sweep time of the television camera. The ver- 
tical deflections of each line represent the tones of grey as 
scanned. Then a ramp signal was added to the horizontal sweep 
lines to separate them out vertically. This ramp signal was 
synched to the vertical sweep time of the television camera. 

Thus a quasi-television picture is made with vertical voltage de- 
flections replacing shades of grey. 

The second method used to display the television scan was to 
select out one individual horizontal scan line with the system as 
shown in Figure 8. The '7>|clock delayed by B clock 1 ' mode of the 
Techtronix 5^5 oscilloscope was used. The A delayed by B mode is 
a system where first a ramp signal set’ to the B clock time base 
is generated. This ramp signal is compared to a value dialed in 
by the scope operator. When the two signals are equal, the A clock 
is armed to take the next trigger it sees and, when triggered, to 
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Figure ?. 
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Figure 8. 





sweep once on its time base 



Thus one sweep (A) in a sequence (3) 



can be selected. By having the B clock synched to the vertical 
svjeep of the television camera and displaying with the A clock . . 
synched to the horizontal sweep of the television camera, a select- 
ed horizontal sweep of the television camera can be displayed. 
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CH APTEi-1 3 

Chapter 3 describes the experimental ly observed results and 
the descriptions of the phase processes inferred from this data. 
3-1 Beam Description 

The individual intensities of the local oscillator beam and 
the transmitted beam after propagation through the atmosphere and 
receiver optics are shown in Figures 10 and 11 respectively. 

These photographs were made by simply blanking off the other beam, 
and so represent typical intensity distributions at the film sur- 
face for each of them. The picture of the intensity over the 
receiver aperture shows the effects of scintillation, ie., a non- 
uniform distribution of intensity with a spatially random pattern. 
The blank spot in the middle and the strip running off to the 
right in Figure 11 are due to the secondary mirror of the Newton- 
ian telescope used as the receiver. Although this mirror blanks 
off an appreciable portion of the photograph, enough of the aper- 
ture is seen to allow analysis. 

When the two beams are present simultaneously and allowed 
to Interfere on the film surface, the fringe patterns shown in 
Figures 12, 13 and 14 are created. It was found through experi- 
mentation that easily discernable changes in the fringe spacing 
over the aperture could not be seen unless the number of fringes 
was reduced to approximately twenty or less. If the model of Ap- 
pendix A is assumed, then equation (A-12) describes the fringe 
pattern as 





Figure 10. Local Oscillator Beam Intensity at Detector 
Surface 




Figure 11. Transmitted Beam Intensity at Detector Surface 
after traversing Atmospheric Path and Receiving 
Optics . 
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where 




( A- 12') 



U x represents the amplitude of the envelope of the field 
sent through the atmosphere 

U 2 represents the amplitude of -the envelope of the local 
oscillator field where a is the phase angle between E, 
2 and E,> * “* x 

N represents the average collected background radiation 
intensity 

L represents the losses in sending a beam through the 
atmospheric channel and in collection 

P represents the direction cosine between the two waves 
p at the detector surface (see Fig. 1) 



y[ represents the impedence of the atmosphere 



U represents the magnif ication of the telescope and eye- 
piece 

e x(P) 

represents the scintillation process due to turbulence 
(jl(p) represents the phase distortion process due to turbulence 
This intensity expression reflects the fact that this pattern is 
assumed to be a time sample of the fields. It contains only spa- 



tially varying terms and implicitly assumes that this pattern is 
taken for a particular time t. At some later time a photograph 



would have another time evoltion of the sample functions X(p) 
and (p( p ) . 

As shown in Appendix B, the geometry can be selected so that 



the fringes occur along one of the axes in the p plane. With this 
geometry the direction cosine, B , can be related to the difference 
in approach angles , H , of the two beams so that the intensity be- 
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Figure 12. Sequence A Fringe Patterns. These consecut 
photographs will be analyzed in the followin 
figures! Figs. 15 

18 

22 

25 
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Figure 13. Sequence B Fringe 
Patterns. These three con- 
secutive photographs will be 
analyzed in the following 
figures! 

Figs. 16 
19 
23 
25 








Figure 14. Sequence C Fringe Patterns. These consecutive 
photographs will be analyzed in the following 
figures! Figs. 17 

20 

24 

27 






comes 
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+ 2uL tJ^ 



( 12 ) 



where 



and 



a = the phase angle between E^ and E ? 
sin£ = |B p j /\S\ c? S. 

£C< 1 



It can be seen that the fringes are caused by the phase modulated 
sinusoid riding on the intensities of the individual signals. 



is small compared to the distance over which significant changes 
occur in the values of <j>( p ) , the phase distortion, individual 
fringe lines will not experience much perturbation. Since a 
maximum of twenty lines distributed over the three inches of 
aperture show appreciable perturbations of the fringe lines, a 
minimum significant distance in<^6(p) can be estimated at (3/20) 
inches = k millimeters. This rough estimate of significant di- 
mensions in(j)(p) conforms to data given, in the literature and 
general experience from longer paths. 

The amplitude scintillation in the three sets of fringe pat- 
terns does not appear to be severe, but this is due to a selection 
process for just this quality. 

Another notable feature of these slides is the spatially vary- 



Figs. 12 and 13 are narrower on the left side of the aperture. 
Especially in the photographs of Fig. 13» one can see that the 



Further, it can be seen that if the period of the sinusoid, 



ing fringe pattern across the aperture. The fringe spacing in 
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fringe spacing noc only varies spatially, but also changes radi- 
cally from one photograoh to the next. The photographs in Fig. 14 
show a. rotation of .Lines between photographs. The photographs of 
Fig. 12 were taken just prior to those in Fig. 13 , while those in 
Fig. 14 were taken on another night without any appreciable modi- 
fication of the system having taken place. 

A photograph of the system through a non-atmospheric path 
could not be obtained unless the system as installed was changed 
significantly. The divergence of the wave emerging from the laser 
or laser— collimator when beamed into the receiver telescope was 
such that its focusing point was moved too far out to use the ca- 
mera eyepiece as mounted. This camera eyepiece was needed to match 
curvatures with the LO wave. Other eyepieces required optics in 
thejl'O path. When the divergence was reduced so that the system 
could be used "as is", the aperture was not fully illuminated, and 
very little of the aperture could be seen around the sedondary mir- 
ror. Thus system effects can not be directly taken out. However, 
the availability of consecutive photographs allows analysis of the 
changes in successive photographs. The out leal system is the same 
between photographs except for advancing the film in the camera, 
so that changes in the fringe spacing would indicate atmospheric 
effects. For this reason, the photographs of Fig. 14 were included 
to obtain another consecutive look git the lumped, system-atmospheric 
effects on a night different from that of Figs. 12 and 13 . 

3-1 Full Aoerture Analysi s 

In an effort to obtain a feel for the processes that occur 
over the entire aperture, the slides were projected on a screen 



The output of the television 
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and scanned with a television camera, 
camera was displayed, on an oscilloscope using the first of the 
metnods described in section 2—4, This method, is essentially 
the construction of a television picture but with the shades of 
grey replaced by "vertical deflections on the individual scan lines, 
|This method gives che appearance of a three dimensional map of the 
intensity over the aperture and the results are shown in Figs. 15 , 
16 and 1?. 

The main interest in displaying in this manner is to get 
some feel for 

(1) Phase changes along the fringes 

(2) Phase changes across the fringes 

(3) Scintillation 

(4) Systematic effects 

The deep troughs on the edges of the photographs in Figs. 

15 > 16 and 17 are the locations of black non-reflecting tape put 
on the projection screen to' aid in centering the picture. In e~ 
very photograph, the fringe spacing increases towards the right. 
This is considered a system constant and. as such will be dis- 
counted in further analysis. 

The fringe lines tend to curl towards the lower right-hand 
corner. This curl is stronger in the lower half of the aperture 
and can be seen particularly well in photograph of Fig, 1?. This 
is an example of the phase distortion along the f^ringe lines. A 
slightly different effect is the amplitude modulation along the 
fringe lines caused by phase distortion. 

The scintillation can be seen best on the peak lines as on 
amplitude modulation. Where the modulation produced by the fringes 
is shallow, generally on the right hand sides of the photographs, 





Figure 15 . Sequence A Full Aperture Display of the 
Fringe Pattern of Photograph <fl and*f 2 
by Method of Section 3-2. 
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Figure 16. Sequence B Full Aperture Display of the 
Fringe Pattern of Photography 1 by Method 
of Section 3-2. 




Figure 17. Sequence C Full Aperture Display of the 
Fringe Pattern of Photograph#! by Method 
of Section 3-2 
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scintillation tends to swamp out the detail of the pattern. 

The analysis of the phase process across the fringe grain 
fill be done in detail in the next two sections. It is worth- 
•?nile noting at this ooint that the fringe pattern is generally 
risible and enough definition exists to justify a look in detail. 
^ ~ 3 Phase Ana lysis 



Using equation (B-2) as the model for this fringe line sys- 
;em, the intensity is 

-if. .2 2 2 2 2 2X(p) 

I(p) = 2^||U 2 | + N + U L juj e 

+ 2uL|U 2 ( |oJ o X P oox[^ . ^( x ,y o ). a ]| (12) 

If a narrow strip in the y dimension is examined, narrow enough 
so that the processes do not change with y, the intensity ex- 
pression becomes 



Kx,y 0 ) = 



U, 



2 2 2 
+ u L 

X(x,y ) 



+ N + u L 1 tM e 



2 2X(x,y ) 



° cos[-J~— -^>(x,y o )-a^ 



+ 2ul|u 2 | JuJ e u cos - <l>(x,y^)-a] l f ( 13 ) 

Assuming that the spatial sinusoid is visible despite the random 
amplitude modulation from scintillation, the distance between suc- 
cessive peaks or valleys can be measured. This indicates the sin- 
usoid argument going through 2tt radians or one period. Thus, 
choosing the peak lines as the point to measure period length on, 
we know over one period 

-cj)( x ,y o )] - [ ^'-<j>(x',y 0 )] = 2„ (14) 

where 



x and x* are the points where two adjacent maximum values 
of the sinusoid are located. In the remainder 
of this thesis x and x* will indicate adjacent 
maxima. 
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Since x and x* can be measured, we get 



<j>(x,y 0 ) -<jxx\.y 0 ) = 2n[~( x-x ' ) - 1] (15) 

rhus over the distance (x-x') we get a measure of how the phase 
distortion process acted. If we continue this analytic procedure 
across the fringe pattern, vre have some notion of how <jpK x ’ y o^ 
lhanges with x. 

It has been shown by Fried that geometrical interpretations 
ian be placed on the form of <p(x,y o ). If a Taylor type expan- 
sion is made of 4>< x,y Q ), we get 



<^>(x s y 0 ) = A + Bx + Cx 2 + . . . (16) 

hitting this expression for (j)(x,y 0 ) into Eq. 13 gives 



I(x » y 0 ) ~ 2 



U, 



2 2 2 
+ N + u L |U 



2 2X(x,yJ 



+ 2uL 



U, uJe 



X(x,y ) 



cos[ 



2tt^ 



a - A - Bx - Cx . . 

(17) 



Jere it can be seen readily that 

A represents a phase origin relative to the position 
coordinates x 

Bx represents a linear phase change, a tilt relative to 
the direction of travel 

2 

Cx ... represents a spherical phase change, representing 
a spherical wavefront on the received wave 

md the coefficients give a good measure of the importance of each 

;erm; the number of terms needed to adequately represent the phase 

>rocess indicate in some sense how simple the phase changes are. 

There are two problems Involved with applying this analysis 

irectly to the fringe patterns made for this thesis. The first 

.s that no non-atmospherlc pattern exists to take out the purely 
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systematic effects. The only method of removing these systematic 
affects is to examine consecutive patterns and remove the constant 
affects; this leads to the second problem. The only distinctive 
feature that can be used for locating the same narrow strip in 
consecutive patterns is the secondary mirror. Thus without taking 
iata for several entire slides, which would take a prohibitive 
amount of time, only limited strips can be handled in each pattern 
for comparison with the next. Fortunately, the regions just above 
md below the secondary mirror are not unique compared to the 
remainder of the fringe patterns in the rest of the aperture, as 
jan be seen by examining Figs. 12, 13 and 14. Thus these regions 
ire a reasonable choice for analysis. 

} - 4 Nar r ow St ri p Ana l ysis 

To examine how the fringe patterns made for this thesis 
changed in phase, narrow strips above and. below the secondary mlr-- 
'•or were measured In detail. The slides containing the fringe pat- 
terns were projected on a screen, then scanned using the televi- 
sion camera. The fringes vie re aligned so that the horizontal 
:can of the television camera was across the grain of the fringe 
pattern. Selected horizontal lines from the television camera 
jrere displayed and the locations of the peaks of the fringe lines 
-ere measured for five to fifteen successive horizontal lines. 

'his corresponds to a strip one to three millimeters wrae over 
.he objective aperture. Typical displays of these horizontal lines 
.re shown in Figs. 18, 19 and 20, where the top line corresponds 
a scan line in the top half of the fringe pattern; the middle 
ine to a scan line in the middle, and the bottom line to a scan 



! 



- 
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Figure 18. Sequence A 
Typical Horizontal Lines 
from Television Camera 
for Photograph 1 



Figure 19. Sequence B 
Typical Horizontal Lines 
from Television Camera 
for Photograph 1 




Figure 20. Sequence C 
Typical Horizontal Lines 
from Television Camera 
for Photograph St 1 
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^ine in the bottom half of the fringe pattern. The relatively 
small modulation in the middle of the middle lines is due to the 
secondary mirror. The relatively deeo excursion of the scan lines 
at the edges is caused by black non-reflecting tape placed on the 
screen co provide a means of centering the fringe pattern on the 
television monitor. The spatial scale of the oscilloscope display 
can be related directly to the spatial scale of the aperture. The 
conversion is 0.40 inches of aperture is equal to 1 centimeter of 
oscilloscope display. This conversion will be applied to all fu- 
ture figures that deal with the aperture. 

Typical unaveraged strips are shown in Fig. 21, where the 
peak lines are drawn in. The locations of these lines are simply 
averaged along the fringes. The distance between the averaged lo- 
cations of the peak lines across the aperture is plotted in Figs. 
22, 23 and 24. From Sq. 15 this distance between peak lines can be 
related to the phase distortion and fringe period as 

(x-x*) = [<j>(x f y 0 ) - 4>(x ' ,y Q ) + 2rr] (18) 

Thus measuring (x-x*) gives a number proportional to (jp over 
the distance (x-x*). Further geometric interpretations about the 
form of (j) can be made as successive (x-x*) values are measured and 
plotted across the aperture. A direct measure of the spatial 
frequency — is not available, however this term can be assumed 
constant across the aperture. Then any changes in the (x-x') 
across the aperture can be attributed to changes in <j? (x,y Q ) and 
any changes in (x-x*) between photographs can also be attributed 
to the change in (|)(x,y o ) with time. 

This method of plotting the average distance between peaks 



ft 30 /C me secor^o^/LY M/2eoZ. 
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Figure 21. Typical Unaveraged. Locations of the Peaks of a 
Fringe Pattern. This example comes from Photo- 
graph 1 of Sequence B. 




over Objective Aperture Sequence C 
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gives a visual display of the general type curvature of the phase 

process across the aperture and seems to snow a general quadratic 

form indicative of spherical distortion of the incoming wave. 

Another method of plotting is to use the averaged location of 

the peak lines on the abcissa and to plot the 2 tt radian change of 

the sinusoid argument between successive peaks on the ordinate. 

This is shown in Figs. 25» 26 and 27. These graphs have a linear 

mean square estimate line plotted for reference. The first aspect 

to note is the dominance of the linear term in every plot. The 

lines have some higher order deformations, but the linear term is 

obviously most important. The second aspect to note is the change 

in slope between photographs in Fig. 26. This means a significant 

change in tilt of the incoming wave has taken place. We can 

equate the argument of the sinusoid in Eq. 13 to the estimate lines 

in Figs. 25s 26 and 2? as 
2rrji c x 

r-y— -<f)(x,y 0 )-a] - k-jX + k g (19) 

where 

\r - XX T1 the linear estimate line slope measured for 
K 1 2 each photograph 
°x 

U,, the covariance between the averaged location of the 
peak lines, x, and the total phase change, y, at 
that x 

2 

o x the variance of the averaged location of the peak lines 

m the mean of the averaged locations of the peak lines 

m the mean of the total chase change across the anerture 

y U 11 

kp = m - m the linear estimate intercept for a 

y photograph 

A 

The change of this argument between successive photographs is 
<£(x,y o ) “ (jy\x\ y Q ) = (k 1 -k , 1 )x - (k 2 ~k’ ? ) 



( 20 ) 
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Figure 27, Total Phase Change Across the Receiving 
Aperture Sequence C 
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where the primes refer to the succeeding photograph. The (k^-k 8 ,,) 
is not meaningful since many steps in the analysis Introduce arbi- 
trary references in phase. However, the change in slope is easily 
measured. These changes in k^ are shown in Table 1 and the equiva 
lent change in tilt of the incoming wave is also listed. 




radians/inch) 



Figure 



Figure 



25 

26 



Above secondary mirror 
Below secondary mirror 
Above secondary mirror 
Below secondary mirror 
Above secondary mirror 



1 . 2r? 

2 rr 

4 . 25 tt/. 75 tt 

8 . 75 n 

.15tt 



Tilt change (in 
seconds of arc) 

3.03 

5.15 

10.9/1.925 

22.5 

.386 



Figure 2? 

Below secondary mirror .75 n 1.925 

Table 1 Measured Changes in Tilt in Three Gets of Consecutive 
Photographs (from Figs. 12, 13 and 14) 

As is discussed in Appendix A, tilt changes less than 15 
microradians ( 3.2 seconds of arc) will be indistinguishable from 
possible mirror tilt caused by vibrations. This limit is exceeded 
only in Fig. 26 . If the photographs used to generate Fig. 26 are 
reexamined (see Fig. 13), it can be clearly seen that the fringe 
spacing changes dramatically from the first to second photogi-aph, 
an effect not seen in Figs, 12 and 14. 

•In an effort to determine the higher order deformations, a 
log-log plot was constructed using the average location of peak 
lines as the abscissa and the phase change s,with the linear esti- 
mate subtracted out, as the ordinate. However, only systematic 
effects are visible. 



CHAPTER 4 



This chapter summarizes the results presented in Chaper 3 , 
and draws some conclusions about the general procedure. Some sug- 
gestions for improvement of this system are also made. 

4 - _ JL C onclusion s 

Using the methods outlined in Chapter 3» it is clearly pos- 
sible to measure phase in a nari'ow strip across the objective 
aperture. Using the resulting phase plots, it is possible to com- 
pute a change in tilt between successive photographs. In this ex- 
periment, the linear terms and systematic effects masked possible 
higher order deformations caused by atmospheric turbulence, but 
with more care on a longer path these effects should be measurable. 
The details of the analysis show that the minimum tilt changes, 
separate from mirror induced tilt, measurable by the present sys- 
tem was approximately 5 seconds of arc. It is interesting to note 
that the measured tilt ranged as high as 22,5 seconds of arc, a 
large value for a three inch aperture and a short path. 

When the fringe pattern produced by the two beams is looked 
at through a Wratten filter, the system dynamics can be seen graph- 
ically. The fringe lines wave r along their length and, a^s their 
spacing gets wider, a bowing travelling across the pattern is of- 
ten seen, A third effect is a shimmering of the entire pattern, 
possibly indicating a time change of phase origin suggested in 
section 3 - 3 » 

The system of oscilloscope display and visual reading is 
clearly a tedious method of data taking. The amount of data 
that can be taken in this way is small, and is not compatible 
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wit h machine computation. Judgement can be exercised to pick 
out the peak lines from scintillation when this display method is 
used * but the peaks and valleys of the fringe pattern are the only 
parts that can be used for visual data taking. 

The lack of a system-only series of photographs was a hin- 
drance to the analysis. Some further feel for the motions of the 
down-range reflector would have been invaluable in interpretating 
the results. 

4-2 Suggestions for Impr ovement 



The first suggestion would be to set up a path in the la- 
boratory long enough so that the divergence of the wave entering 
the telescope is not too fast. At the same time, the beam should 
illuminate a large amount of the aperture. Thus, using the same 
system as one would use when taking atmospheric data, a recording 
of the system alone could be made. This would be most helpful in 
analyzing single photographs, and sequences would not be necessary. 

A second suggestion is to convert the video output of the 

television camera to digital data for later machine analysis. The 

basic equipment is available, and only a gating circuit and sam- 
logic 

pli*)g sequence/need be added. The master clock that generates 
the pulse train used to synchronize the television camera runs at 
110,250 Hz; the first frequency used in the television synchroni- 
zation is at 1 575° Hz. The sampling pulses could be selected to 
run off this master clock with some flexibility as to the desired 
frequency, and synchronization with the television camera could 
still be retained. With the data in digital form on tape, a di- 
gital computer could be used to analyze the data, allowing a de~ 
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tailed analysis over the entire aperture. 

As longer paths are used, the problem of reflector induced 
tilt would be lessened. The increase in beam diameter would not 
be too large as path length is extended, but beam motion due to 
reflector tilts would be more apparent. The motion of the beam 
could be observed and data could be taken only when the beam is 
steady on the target. Thus mirror induced tilt would be minor. 
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APPENDI X A 

This appendix contains an analysis of the heterodyning sys- 
tem in more detail than in Chapter 2, a fuller explanation of the 
data analysis system, and electronic circuit diagrams, 

A - 1 Heterodyne S y stem Analysis 

The scheme outlined in Figure 2 for creating a spatial in- 
terference pattern proportional to the phase distortion process 
caused by atmospheric turbulence can be analyzed using simple 
models . 

The front end of the laser (typically the output of the 
Spectra-Physlcs Model 119 laser would be 0,1 milliwatt at the 
front end and 10 microwatts at the back end) is collimated and 
transmitted through the turbulent atmosphere and then reflected 
back over the same path by the down range reflector. After col- 
lection by the receiving telescope this signal is mixed with a 
local oscillator signal taken from the back end of the laser. 

These two waves can be approximately described before transmission 
as plane waves of finite spatial extent. As the finite spatial 
extent of the fields should not cause important effects for this 
experiment, it is not explicitly included. Thus in complex nota- 
tion 




-JP-L* r jwt 



(A-l) 



e 





r jwt 



e 



where 




\ = 6328ft 



(A-2) 
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The atmospheric effects on the propagating beam include dif- 
fraction, scattering, and absorption losses; atmospheric turbulence 
modulation; beam spreading due to imperfect reflector-mirror flat- 
ness; phase modulations due to mirror motions caused by vibra- 
tion, Collection of the returned beam will introduce further 
losses due to the finite collecting aperture which is smaller than 
the beam size. Background radiation will also be collected in the 
receiving telescope. The various effects will be described as fol- 
lows* _ 

Y(r,t) 

e represents the turbulent channel 



where 

Y(r,t) = X(r,t) + j<j;>(r,t) 

X(r,t) represents the amplitude scintillation process 
4>(r,t) represents the phase distortion process 

L(r,t) represents losses due to 

1. Diffraction 

2. Particle scattering 

3. Molecular absorption 

4. Mirror imperfections 

_ 5. Collection losses 

jQ ^gp^sents phase modulations due to mirror movements 

E (r,t)e n represents the wide bandwidth Gaussian noise 

coilected by the receiver 

Using the channel model of Chapter 1, the transmitted field will be 



received as 



£R e c (r > t) 



_ jfVr Juot Y(r,t) M(r, 

L(r,t)U.,e e_ e e 

„ J&(r,t) Jwt 

+ E (r,t) e e 



t) 



(A-3) 



The collecting telescope will minify the portion of this re- 



ceived beam within its aperture. As long as one looks "close 
enough" to the eyepiece, the minified, image is a clean represen- 
tation of illumination in the object space. (Going further away 
gets into the Fresnel zone of the secondary mirror and aperture 
edge effects.) The collected field is mixed with a local oscil- 
lator field by a beam splitter and projected on a plane screen for 
detection. The geometry of this plane screen (Fig. 1) is selected 
so that z=0 on the surface; the transmitted wave approaches the 
surface normally and the local oscillator signal approaches the 
surface obliquely, and thus has a direction cosine on the surface. 
The fields using the geometry of Figure 1 are 

„ -j'P'p’ jwt 

E Lo (P*t)=U 2 e e (A-4) 

Y ( 0 • t ) iuut 

-Trans (p$t) = L(p,t)uU 1 e e 

j© (P,t) jwt 

+ pS n (p,t)e e (A-5) 

-total ~ ~Lo^ p?t ^ * -Trails^’ ^ (A~6) 

where 

p* = the two-dimensional vector in the plane of the 
detector surface 

M = the collecting telescope-eyepiece raagnif ication 

| 

The detector, in this case film, will measure the intensity which 
T , v ^ ^ total l [ ^t otal ^ 

Further, a black box model of this detection process includes a 
time averaging over the period of exposure as shown in Fig A-l. 

T is the period of the shutter opening and is on the order of 2-10 
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1 second* 



i»c ss is 



for the 35« camera. The 1M» 
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u 



\(p,t) 



2 X(p,t) 

2L(p,t)u U^E n ( p , t ) c cos[© n (p,t) - 4>(p,t)- tt(p,t)] 



and 2aU 2 E n (p» t)cos[p • p -© n (p',t)] 

The integral of the first term, the square of the envelope 
of a complex Gaussian process, is a relatively long-term integra- 
tion of a very wide bandwidth process; thus the time integration 
can be replaced with the ensemble average 

* tf 

p 

where N' is proportional to the average power in the relevant 
background radiation. 

Similarly, the time integrals of the terms containing both 
E n ( p , t ) and @ ( p,t) may be replaced with ensemble averages since 
both E n (p,t) and0^p,t) are very wide bandwidth processes com- 
pared to the integration time T. Thus, for example 

~ | ZJLi (7/ E n <f,t) Cos Zfe f- Cf.ojjt = f 

Further, the two processes E n (p,t) and © n (p,t) are independent, 
so that 

1 1 ©»<?,<£)] = £[**■[¥! 



f" >•/'- cz,sjj e f~ e>n<f,t>] J 



Also,0 n (p,t) is reasonably modeled as approximately uniformly 
distributed between 0 and 2rr at any time during the period 0 to T. 
Thus E £cos[B p *p ~0 n (p,t)]^ tends tovrards zero and the time inte- 
gral is anproximately zero over the entire surface. Similarly, the 
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term tends towards zero. 

This integration model of the detector implies that the de- 
tector records energy. For this model, since the processes are 
slow compared to the shutter time, a typical intensity, I("p), has 
been used with the understanding that multiplication by T, the 
integration time constant, will yield the proper energy terms. 

The resulting output from this integration model is thus 



The loss term L(p) is reasonably modeled as a uniform quan- 
tity over the strong portion of the transmitted beam since scat- 
tering and molecular absorption are essentially non-pref erential 
and are less apparent in the strong portion of the beam. As it 
is this portion of the beam that will be collected, the term L( p ) 
will be replaced with L. 

The term ty(p) represents the phase modulation due to mirror 
vibrations further acted on by the turbulent channel. Due to 
mirror mount massiveness, the value of ty(p,t) is constant with 
time over the shuttering period. The spatial variation of ty(p) 
is critical, however, since the mirror movements that can occur 
cause a tilt on the reflected wave. This tilt is of the same form 
that one expects as the first order effect of atmosphere dis- 
tortion. 

There are arguments that seem to Indicate that while this 

tilt caused by the mirror is a limitation on the system, the mea- 




sureraents of the atmospheric effects can still be made. The first 
is that tyoical tilt changes between photographs are on the order 

■ 

; of 15 microrad ians . If this tilt change were due only to mirror 
i movements, this would mean a beam movement of 0. 05 inch at the 
telescope. The overlap of the returned beam over the receiving 
was approximately 0. 5 inches, meaning that this movement could 
not be seen. However, more significant tilt changes would be 100 
plus microradians, which translates to a 0.5 inch travel at the 
telescope. Thus this movement would be seen and neglected. Thus 
it is recognized that mirror motions are critical, but for this 
path do not swamp out larger values of atmospheric tilt. There- 
fore, ty(p) will be lumped in 
as 



I(P) = 



Jl 

?S\ 



N 



(j)(p), and the intensity can be written 
2 2 



-!- N + u. L |U 1 | 
X(p) 



2 2 X(p) 
e 



2 uL N U^je cos C 8 p . P -^>(p)-tt]j (A- 12 ) 



where the intensity contains a terra phase modulated by the atmos- 
pheric phase distortion. 

A - 2 Electronics 



The data analysis system was built around the Packard Bell 
920 television camera and is shown in block diagram form in Fig* 
A- 2 . As produced, the camera has a 2*1 random interlace which 
was altered to a single vertical scan to allow accurate reading 
of the fringe pattern. The sweep of the television camera was 
made into 256 horizontal lines each vertical sweep, with synchron- 
ization between the horizontal and vertical rasters. The circuits 
used for this system were built around RTL digital integrated- 



circuits 
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Figure h-% Data Analysis System Block Diagram 
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A master clock shown in Fig. A-*3 was used to generate the 
pulse train for the system. 

A series of counters and pulse shapers were used to generate 
synchronization pulse trains from the master clock signal. In 
each case a synchronous counter was used to minimize countdown 
delays inherent in sequential counters. Figs, A-4, A- 5 and A-6 
are circuit diagrams for these counters and pulse shapers. In 
each of the circuits the counting is done by the JK flip-flops 
gated on and off by the logic. The shapers were pulse generators 
with a variable resistor used to set the output pulse width. 

The synchronization pulses were put into Schmidt trigger 
circuits to provide some power gain and also to buffer the count- 
er circuits. 

Using the system outlined above, the television camera scan 
was made to cover the 2$6 horizontal lines in each vertical sweep 
with no interlace. The starting point for each vertical sweep 
was the .same so that the video output could be displayed on an 
oscilloscope, using the synchronization pulses as triggers, with 
precision and stability. 

Before display, the video output of the television camera 
was put into a low pass filter with a 3 db point at 4 MHz. The 

low pass filtering was done to remove a large amount of high 

* 

frequency hash. Fig. A-? shows the improvement of the display pro- 
duced by this filtering. 






* 
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Figure A— ^ Master ClocM and Pulse Generator for 
Synchronization System 




Figure A™^. Eight Counter and Pulse Generator 
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Figure A~& Seven Counter and Pulse Generator 





Figure A-$ Four Counter and Pulse Generator 
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Figure A-7. Improvement in Picture Quality from Insertion 
of Low Pass Filter at the Television Video 
Output. Upper Picture is without filtering; 
Lower Picture is with filtering. 
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APPENDIX B 



This appendix contains a discussion of the data analysis sys- 
tem limitations. 

B - 1 Fringe Pattern Model 



As shown in Appendix A, the intensity for plane wave models 



is 

_ n r 2 2 22 2 2x( p ) 

1(p > = 41 '- 2l +N + ^ ) l- u -il e 

+ 2ULJUJ | U 2 ] e oos[B p - P - tj>( pl-aH (A-12) 

The fringes result from the spatial sinusoid, which is proportional 
to the product of the two signals. The fringes ride a bias which 
is the sum of the intensities of the two signals. The terms re- 
presenting atmospheric turbulence effects can be thought of as 
minified representations of the processes at the objective of the 
receiving telescope. Thus 

x(p) = x(p'/u) 

4>Cp ) = <$>(?• /u) 

where 



p is the tvro dimensional vector at the detector surface 

p’ is the two dimensional vector at the objective plane 
oriented so that the geometry, when propagated through 
the receiver optics, matches that of p. 

U represents the magnification of the telescope-eyepiece 
combination 

Further, aligning the geometry of this model so that the fringe 
lines coincide with one of the axes in the p plane simplifies the 
notation. With this assumption, the intensity is 



2 
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Kp) 




2 22 2 2X( p ) 

+ N + U L luj e 



+ 2uL^U 1 | j U 2 |cos[B x x -<j>(p)-a] 

where 

8 x = | B |^x/l 8 l ) = |B[sinjJ,= 2n l/\ 

t represents the angle between the two interfering waves 
and is set by the experimenter so that £<6<1 

By looking at a very narrow strip in the y dimension so that vari- 
ations in y do not have a chance to occur, the intensity model can 
be written as 



I(P) 




22 2 2X(x v y ) 

+ M L |.U 2 | e 



+ ZuI^Uj lu^ X ’ yo) oos[~Jk - (j>(x,y 0 )-tt]| ( B-l ) 

The problem of seeing the fringe pattern is related to the 
depth of modulation, which can be shown in this model to depend on 
the visibility ratio V, defined as 



Y - amplitude of sin us oid 
amplitude of bias"* 

, . . , X(p) 

= 2uL K 1 Mi 

. . 2 " 2 2 2. ,2 2X(p) (B-2) 

U + N + u £ |UJ e 

For this experiment, N was x negligible compared to the other terms 
in the denominator and will be neglected. The other terms in the 
denominator represent the local oscillator and received beam in- 
tensities. .The shadow pattern of these intensities can be seen 
during experimentation. From observation, the shadow pattern of 
the received field is brighter than the local oscillator by an 
estimated factor of three to four. On a longer path, losses will 
surely be more severe and this ratio would be reversed. If Eq. B-2 
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is manipulated slightly, and the inforraation on the ratio of 



the two beam intensities incorporated, the visibility ratio be- 



comes 




= 2 



1 + K Z 7^ 



1 f 4e 2X(l ’ T 



(B-3) 



It can be seen that the scintillation process is the fundamental 
limitation on modulation depth for this path. This is plotted in 
Fig. B-l. 




Figure B-l Plot of Visibility Ratio with changes in 
the value of scintillation 



A further limitation peculiar to this particular data 
analysis system was the resolution available using an oscillo- 
scope display for detailed analysis. The scope face was ten 
centimeters across and distinctive portions of the waveform 
could be located accurately to one-half millimeter as long as 
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the surrounding clutter was not too bad. However, at least a 
three to five millimeter separation was required before it was 
possible to reliably distinguish peaks from clutter for the typi- 
cal depths of modulation produced by the scanned fringe lines. 
This means that 20 to 33 lines maximum could be handled on the 
scope. It was found that (see section 3~i) the effects of the 
atmospheric turbulence on the fringe perturbations did not become 
easily visible until the number of fringes across the aperture 
became twenty or less. Thus the use of the scope was not a 
serious limitation for this experiment. 
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